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ABSTRACT 



Context. We report the detection by the AGILE satellite of a rapid y-ray flare from the source 1AGL J151 1-0908, associated with the powerful 
y-ray quasar PKS 1510-089, during a pointing centered on the Galactic Center region from 1 March to 30 March 2008. This source has been con- 
, tinuosly monitored in the radio-to-optical bands by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT). 
OO ■ Moreover, the y-ray flaring episode triggered three ToO observations by the Swift satellite in three consecutive days, starting from 20 March 
"^J" ■ 2008. 

C*") ' Aims. The quasi-simultaneous radio-to-optical, UV, X-ray and y-ray coverage allows us to study in detail the multifrequency time evolution, the 
, spectral energy distribution of this source and its theoretical interpretation based on the synchrotron and inverse Compton (IC) emission mecha- 
. nisms. 

' Methods. During the radio-to-optical monitoring provided by the GASP-WEBT, AGILE observed the source with its two co-aligned imagers, the 
^ j Gamma-Ray Imaging Detector (GRID) and the hard X-ray imager (SuperAGILE), sensitive in the 30 MeV-30 GeV and 18-60 keV energy bands, 
. . respectively. 

^ ■ Results. In the period 1-16 March 2008, AGILE detected y-ray emission from PKS 15 10-089 at a significance level of 6.2-<x with an average flux 
k>( ' over tlle entrre period of (84 ± 17) x 10~ s photons cnT 2 s _I for photon energies above 100 MeV. After a predefined satellite re-pointing, between 
rS | 17 and 21 March 2008, AGILE detected the source at a significance level of 7.3-tr, with an average flux (E > 100 MeV) of (134 + 29) x 10~ 8 
, photons cirr 2 s -1 and a peak level of (281 ± 68) x 10~ 8 photons cirr 2 s -1 with daily integration. During the observing period January-April 2008, 
. . . i the source also showed an intense and variable optical activity, with several flaring episodes and a significant increase of the flux was observed at 
millimetric frequencies. Moreover, in the X-ray band the SvW/f/XRT observations seem to show an harder- when-brighter behaviour of the source 
spectrum. 

Conclusions. The flat spectrum radio quasar PKS 1510-089 showed strong activity between January and April 2008, with episodes of rapid 
variability from radio to y-ray energy bands, in particular with a rapid y-ray flaring episode. The spectral energy distribution of mid-March 2008 is 
modelled with a homogeneous one-zone synchrotron self Compton (SSC) emission plus contributions from inverse Compton scattering of external 
photons from both the accretion disc and the broad line region. Indeed, some features in the optical-UV spectrum seem to indicate the presence of 
Seyfert-like components, such as the little blue bump and the big blue bump. 

Key words, gamma-rays: observations - mechanism: non-thermal - quasars: individual (PKS 1510-089) 



1. Introduction 



Send offprint requests to: F. D'Ammando, e-mail: 
f ilippo . dammando@iasf-roma . inaf . it 

* The radio-to-optical data presented in this paper are stored in the 
GASP-WEBT archive; for questions regarding their availability, please 
contact the WEBT President Massimo Villata. 



Blazars are a subclass of AGN characterized by the emission 
of strong non-thermal radiation across the entire electromag- 
netic spectrum, from radio to very high energies. Their obser- 
vational properties include irregular, rapid and often very large 
variability, apparent super-luminal motion, flat radio spectrum, 
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high and variable polarization at radio and optical frequencies. 
These features are interpreted as the result of the emission of 
electromagnetic radiation from a relativistic jet that is viewed 
closely aligned to the line of sight (Blandford & Rees 1978, Urry 
& Padovani 1995). The Spectral Energy Distribution (SED) of 
blazars is typically double-humped with a first peak occurring 
in the IR/optical band in the so-called red blazars (including 
Flat Spectrum Radio Quasars, FSRQs, and Low-energy peaked 
BL Lacs, LBLs) and at UV/X-rays in the so-called blue blazars 
(including High-energy peaked BL Lacs, HBLs). The first peak 
is commonly interpreted as synchrotron radiation from high- 
energy electrons in a relativistic jet. The second component of 
the SED, peaking at MeV-GeV energies in red blazars and at 
TeV energies in blue blazars, is commonly interpreted as in- 
verse Compton (IC) scattering of seed photons by relativistic 
electrons (Ulrich et al., 1997), although a different origin of high 
energy emission has been proposed in hadronic models (see e.g., 
Bottcher 2007, for a recent review). 

PKS 1510-089 is a nearby (z=0.361) radio-loud highly po- 
larized quasar (HPQ) already detected also in the MeV-GeV en- 
ergy band by the EGRET instrument on board CGRO (Hartman 
et al. 1992). The broadband spectrum of the source is representa- 
tive of the class of FSRQ with the inverse Compton component 
dominated by the y-ray emission and the synchrotron emission 
peaked around IR frequencies, even if it is clearly visible in this 
source a pronunced UV bump possibly due to the thermal emis- 
sion from the accretion disk (Malkan & Moore 1986, Pian & 
Treves 1993). 

PKS 15 10-089 has been extensively observed and studied in 
the X-ray band by the satellites EXOSAT (Singh, Rao & Vahia 
1990, Sambruna et al. 1994), GINGA (Lawson & Turner 1997), 
ROSAT (Siebert et al. 1998), ASCA (Singh, Shrader & George 
1997) and Chandra (Gambill et al. 2003). The observed X-ray 
spectrum was very flat in the 2-10 keV band with photon index 
of T ^ 1.3, but it was steeper (r 1.9) in the ROSAT band- 
pass (0.1-2.4 keV), suggesting the presence of a possible spec- 
tral break around 1-2 keV. The difference of the photon indices 
could be due to the presence of a soft X-ray excess. Observations 
by BeppoSAX (Tavecchio et al. 2000) confirm a possible pres- 
ence of a soft X-ray excess below 1 keV. Evidences of a similar 
soft X-ray excess has been detected in other blazars such as 3C 
273, 3C 279, AO 0235+164 and 3C 454.3, even if the origin of 
this excess is still an open issue. 

During August 2006, PKS 15 10-089 was observed in a rela- 
tively bright state by S uzaku over approximately three days and 
with Swift 10 times during 18 days as a Target of Opportunity 
(ToO) with a total duration of 24.3 ks. S uzaku measured a very 
hard X-ray spectrum (r < 1 .5), which seems to exclude models 
in which X-rays are produced by synchrotron radiation of the 
secondary ultrarelativistic population of electrons and positrons, 
as prediced by hadronic models. The SwiftPiRT observations in- 
stead revealed significant spectral evolution in the 0.3-10 keV 
energy band on timescales of one week, with the spectrum that 
becomes harder as the source gets brighter (Kataoka et al. 2008). 

Gamma-ray emission from PKS 1510-089 was detected 
several times by EGRET with a integrated flux above 100 MeV 
between (13 ± 5) and (49 ± 18) x 10~ 8 photons crrT 2 s _1 and an 
energy spectrum, integrated over all the EGRET observations, 
modelled with a power law with photon index F = 2.47 ± 0.21 
(Hartman et al. 1999). 

In August 2007, AGILE detected an intense y-ray activ- 
ity from PKS 1510-089. In particular, during the period 28 
August - 1 September 2007 the average y-ray flux observed was 
F £ >iooM e v = (195 + 30) x 10~ 8 photons crrT 2 s" 1 (Pucella et al. 



2008). Recently, this source was detected again during high y- 
ray activity states by both the Large Area Telescope (LAT) on 
board the Fermi GST (Ciprini et al. 2008, Tramacere 2009 and 
Cutini et al. 2009) and AGILE (D'Ammando et al. 2009, Pucella 
et al. 2009 and Vercellone et al. 2009). The results of the AGILE 
observations during March 2009 will be published in a forth- 
coming paper (D'Ammando et al., in preparation). 

In this paper we present the analysis of the AGILE data ob- 
tained during the observations of PKS 1510-089 from 1 March 
2008 to 30 March 2008 (see also D'Ammando et al. 2008a). 
We also present the radio-to-optical monitoring of the GASP- 
WEBT during the period January-April 2008, and the results of 
three Swift ToO carried out between 20 and 22 March 2008. 
This broadband coverage over the entire electromagnetic spec- 
trum allows us to build and study in detail the spectral energy 
distribution of the source. 

The paper is organized as follows. Section 2 describes 
the AGILE observations, and the corresponding data analysis. 
Section 3 introduces the Swift data and the relative analysis. 
Section 4 is dedicated to the results of the GASP-WEBT ob- 
servations, while in Section 5 we discuss the spectral energy 
distribution, its implication for the emission mechanisms of the 
source, and finally we draw our conclusions. 

Throughout this paper the quoted uncertainties are given at 
the 1-cr level, unless otherwise stated, and the photon indices are 
parameterized as N(E) oc E~ r (ph crrT 2 s _1 keV -1 or MeV -1 ) 
with r = a + 1 (a is the spectral index). We adopt a luminosity 
distance of d L = 1915 Mpc for PKS 1510-089, assuming z = 
0.361 and a ACDM cosmology with H = 71 Km s" 1 Mpc" 1 , 
Q m = 0.27 and Q A = 0.73. 



2. AGILE Data 

2.1. Observation of PKS 1510-089 

The AGILE satellite (Tavani et al. 2008a,b) is an Italian Space 
Agency (ASI) Mission devoted to high-energy astrophysics, 
with four active detectors capable of observing cosmic sources 
simultaneously in X-ray and y-ray energy bands. 

The Gamma-Ray Imaging Detector (GRID) consists of a 
combination of a pair-production Silicon Tracker (ST; Prest 
et al. 2003, Barbiellini et al. 2001), sensitive in the energy 
range 30 Me V-30 GeV, a non-imaging CsI(Tl) Mini-Calorimeter 
(MCAL; Labanti et al. 2009) sensitive in the 0.3-100 MeV en- 
ergy band, and a segmented Anti-Coincidence System (ACS) 
made of plastic scintillator layers which surrounds all active de- 
tectors (Perotti et al. 2006). A co-aligned coded-mask hard X- 
ray imager (SuperAGILE; Costa et al. 2001, Feroci et al. 2007) 
ensures coverage in the 18-60 keV energy band. 

The AGILE observations of PKS 1510-089 were performed 
from 1 March 2008 12:45 UT to 21 March 2008 2:04 UT, for 
a total of 211 hours of effective exposure time. In the first pe- 
riod, between 1 and 16 March, the source was located ~ 50° 
off the AGILE pointing direction. In the second period, between 
17 March and 21 March, after a satellite re-pointing, the source 
was located at ~ 40° off-axis. Finally, after a gap of 4 days of 
observation due to technical maintenance of the satellite, the 
source was observed at ~ 50° off axis between 25 March 13:09 
UT and 30 March 10:29 UT. Unfortunately during the observa- 
tion the source was substantially off-axis in the field of view of 
SuperAGILE. 
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2.2. Data reduction and analysis 

AGILE-GRID data were analyzed, starting from the Level- 
1 data, using the AGILE Standard Analysis Pipeline (see 
Vercellone et al. 2008 for a detailed description of the AGILE 
data reduction). An ad-hoc implementation of the Kalman Filter 
technique is used for track identification and event-direction re- 
construction in detector coordinates and subsequently a quality 
flag is assigned to each GRID event, depending on whether it is 
recognised as a confirmed gamma-ray event, a charged particle 
event, a single-track event, or of uncertain nature. 

After the creation of the event files, the AGILE Scientific 
Analysis Package can be run. Counts, exposure, and Galactic 
background y-ray maps were generated with a bin-size of 0.25°x 
0.25° for photons with energy E > 100 MeV. 

To reduce the particle background contamination, we se- 
lected only events flagged as confirmed y-ray events, and all 
events collected during the South Atlantic Anomaly were re- 
jected. We also rejected all the y-ray events whose reconstructed 
directions form angles with the satellite-Earth vector smaller 
than 80°, in order to reduce the y-ray Earth albedo contami- 
nation. The most recent version (BUILD- 16) of the Calibration 
files at the time of writing ^publicly available at the ASI Science 
Data Center (ASDC) sitqj and of the y-ray diffuse emission 
model (Giuliani et al. 2004) is used. We ran the AGILE max- 
imum likelihood procedure with a radius of analysis of 10°, on 
the whole observing period, in order to obtain the average flux in 
the y-ray band and estimate the diffuse parameters used also for 
measure the daily fluxes, according to the procedure described 
in Mattox et al. (1993). 

2.3. Results 

During the period 2008 March 1-16, AGILE-GRID detected y- 
ray emission from a position consistent with the powerful y-ray 
quasar PKS 1510-089 at a significance level of 6.2-cr with an 
average flux over the entire period of (84 + 17) x 10~ 8 photons 
cm' 2 s _1 for photon energies above 100 MeV. 

Instead, in the period 2008 March 17-21, AGILE detected y- 
ray emission from a position consistent with the source at a sig- 
nificance level of 7 .3-cr. The AGILE 95% maximum likelihood 
contour level baricentre of the source is I = 351.49°, b = 40.07°, 
with a distance between this position and the radio position 
(I = 351.29°,/? = 40.14°) of 0.17°. The overall AGILE error 
circle, taking both statistical and systematic effects into account, 
has a radius r = 0.50°. The average flux above 100 MeV during 
this second period, with the source located ~ 40° off the AGILE 
pointing direction, was (134 + 29) xlO^ 8 photons cm s . The 
peak level of activity with daily integration was (281 +68) x 10~ 8 
photons cirT 2 s _I , showing an increase of a factor two in one day 
and at least three in two days, as the source had not been detected 
for some days after two episodes of medium intensity. After the 
sudden increase, the flux rapidly decreased around March 19, 
2008. 

Fitting the data relative to the period March 17-21 with a 
simple power law model we obtain a photon index of F = 1.81 
+ 0.34. This photon index is calculated with the weighted least 
squares method, considering for the fit three energy bins: 100- 
200 MeV, 200-400 MeV and 400-1000 MeV. The photon index 
obtained for this second period is consistent within the errors 
with the one observed by AGILE in August 2007 (F = 1.98 + 
0.27). 
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Fig. 1. AGILE-GRID y-ray light curve between 1 and 21 March 
2008 at 1-day or 2-day resolution for E > 100 MeV with fluxes in 
units of 10~ 8 photons cirT 2 s _1 . The downward arrows represent 
2-cr upper limits. 



Figure 1 shows the y-ray light-curve between 1 and 21 
March 2008 with 2-day resolution for the first period and 1-day 
for the second period, for photons of energy above 100 MeV. 
The downward arrows represent 2-cr upper limits. Upper limits 
are calculated when the analysis provides a significance of de- 
tection < 3-cr (see Mattox et al. 1996). 

Finally, in the third period between 25 and 30 March 2008 
the source was not detected by the GRID and an upper limit with 
95% confidence level of 54 xl0~ 8 photons cirT 2 s _1 is provided. 

During August-October 2008, Fermi-LKI detected the 
source with an average flux for E > 100 MeV of (55.8 ± 3.3) 
x 1 0~ 8 photons cirT 2 s~ 1 and a peak of intensity of ( 1 65 .9 ± 1 1 .7) 
xlO 8 photons crrT 2 s _1 (Abdo et al. 2009). The peak of y-ray 
emission corresponds to the first flare observed by FerraZ-LAT 
at the end of September 2008 (Tramacere 2008); the average 
flux value confirms the flaring state observed by AGILE in mid- 
March. 

Moreover in August-October 2008, Fermi-LAST observed a 
softer photon index for this source, T — 2.48+0.05 (Abdo et al. 
2009), but this value corresponds to an average value over three 
months of observation in which the source flux was variable, 
whereas the value reported by us refers to a rapid flaring episode. 
The value obtained by Fermi-LAST is very similar to that mea- 
sured by EGRET averaging over all the observations (r = 2.47 
+ 0.21), confirming that the average spectral indexes are softer 
than those measured during short flaring states. The difference 
between the value obtained by AGILE and Fermi could also be 
partially due to the different bandpasses of the two intruments. 



3. SWIFT observations 

The NASA Swift gamma-ray Burst Mission (Gehrels et al. 
2004), performed three ToO observations of PKS 1510-089 in 
three consecutive days with the first occurring on 20 March 
2008. The three observations were performed using all three 
on-board experiments: the X-ray Telescope (XRT; Burrows et 
al. 2005, 0.2-10 keV), the UV and Optical Telescope (UVOT; 
Roming et al. 2005, 170-600 nm) and the coded-mask Burst 
Alert Telescope (BAT; Barthelmy et al. 2005, 15-150 keV). The 
hard X-ray flux of this source is below the sensitivity of the BAT 
instrument for so short exposure and therefore the data from this 
instruments will not be used. 



3.1. Swift/XRT Data 

The XRT data were processed with standard procedures 
(xrtpipeline vO.12.0), the filtering and screening criteria were 
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Table 1. Observation log and fitting results of Swift/XKT observations of PKS 1510-089. Power law model with Nh fixed to 
Galactic absorption. " Unabsorbed flux. b Cash statistic (C-stat) and percentage of Montecarlo realizations that had statistic < 
C-stat, performing 10 4 simulations. 



Observation 
Date 


Exposure Time 
(sec) 


Counts 
(0.2 - 10 keV) 


Flux 0.3-10 keV 

erg cnT 2 s" 1 


Photon index 

r 


X 2 red (d.o.f.)/ C-stat (%)» 


20- Mar-2008 

21- Mar-2008 

22- Mar-2008 


1961 
1966 
1885 


306 
189 
261 


12.20+; " x 10 12 
8.77 + ; | x 10-' 2 
9.48!|fj x 10~ 12 


1.16 + 0.16 
1.53 + 0.17 
1.41+0.19 


0.87 (13) 
467 (48.2)* 
1.34(11) 



applied by means of the FTOOLS in the Heasof t package v6.5. 
Given the low rate of PKS 1510-089 during the three observa- 
tions (< 0.5 count s _1 in the 0.2-10 keV range), we only consid- 
ered photon counting (PC) data for our analysis, and further se- 
lected XRT grades 0-12 (according to Swift nomenclature, see 
Burrows et al. 2005). No pile up correction was necessary. The 
ancillary response files were generated with the task xrtmkarf, 
applying corrections for the PSF losses and CCD defects, and 
we used the latest spectral redistribution matrices (RMF, vOll) 
in the Calibration Database maintained by HEASARC. The 
adopted energy range for spectral fitting is 0.3-10 keV and all 
data were rebinned with a minimum of 20 counts per energy bin 
to use the x 2 minimization fitting technique. An exception is the 
observation of 2008 March 21, when the number of counts was 
so low that the Cash statistic (Cash, 1979) on ungrouped data 
was used. S w//>/XRT uncertainties are given at 90% confidence 
level for one interesting parameter, unless otherwise stated. 

Spectral analysis was performed using the XSPEC fitting 
package 12.4.0 (Arnaud et al. 1996). We fitted the spectra with a 
power law model with Galactic absorption fixed to N H = 6.89 
x 10 20 crrT 2 (Kalberla et al. 2005). Table 1 summarizes the 
most important information on XRT observations and the rel- 
ative spectral fit parameters. 

A variability of about 30% in the X-ray flux of the source 
was observed on time scale of one day. Notwithstanding the un- 
certainties due to the errors on fluxes and photon indexes, the 
XRT data seem to indicate that the X-ray spectrum becomes 
harder when the source gets brighter, confirming the behaviour 
already observed in this source by Kataoka et al. (2008) during 
the S wift/XRT observations carried out in August 2006. This 
is a trend often observed in HBL (see e.g Massaro et al. 2008, 
Tramacere et al. 2007, Kataoka et al. 1999), but quite rare for 
quasar-hosted blazars such as PKS 1510-089. For instance, 3C 
454.3 shows approximately the same spectral slope in different 
brightness states (see e.g. Raiteri et al. 2007; Raiteri et al. 2008). 

3.2. Swift/UVOT Data 

During the three Swift pointings, the UVOT instrument (Poole 
et al. 2008) observed PKS 1510-089 in all its optical (V, B, 
and U) and UV (UVW1, UVM2, and UW2) photometric 
bands. Data were reduced with the uvotmaghist task of the 
HEASOFT package. Source counts were extracted from a cir- 
cular region of 5 arcsec radius, centred on the source, while the 
background was estimated from a surrounding annulus with 8 
and 18 arcsec radii. In the first two days only one exposure per 
filter was available, while three exposures per filter were ac- 
quired in the last day. With the only exception of UVM2, the 
source brightness resulted quite stable in all the UVOT bands, 
with variations of a few hundredth of mag, well inside the typ- 
ical UVOT data uncertainty of 0. 1 mag due to both systematic 
and statistical errors. Average values are: V = 16.94, B = 17.19, 



U = 16.31, UVW1 = 16.64, and UVW2 = 16.55. The UVM2 
frames have low signal-to-noise ratios, thus the source magni- 
tude in this band presents a larger dispersion; the average value 
is UVM2 = 16.47 ±0.14. 



4. Radio-to-optical observations by the GASP 

PKS 1510-089 is one of the 28 y-ray-loud blazars that are 
regularly monitored by the GLAST-AGILE Support Program 
(GASP; Villata et al. 2008) of the Whole Earth Blazar Telescope 
(WEBT). Optical and near-IR data are collected as already cali- 
brated magnitudes, according to a common choice of photomet- 
ric standards from Raiteri et al. (1998). Radio data are provided 
as calibrated flux densities. The reference optical band for the 
GASP is the R band; the corresponding light curve in January- 
April 2008 is shown in the top panel of Fig. 1, with the data pro- 
vided by the following observatories: Abastumani, Calar AltcQ, 
Crimean, Lowell (Perkins), Lulin, Roque de los Muchachos 
(KVA and Liverpool), San Pedro Martir, St. Petersburg, Torino. 
The source showed intense activity during all the considered pe- 
riod, with several episodes of fast variability. At the beginning 
of the optical observing season, the January observations indi- 
cate that the source was in a faint state, around R = 16.6. A fast 
brightening of ~ 1.3 mag in 8 days led the source to R = 15.3 
on February 15. This was followed by a ~ 0.6 mag dimming 
in 4 days. Other minimum brightness states were observed on 
March 23 and in late April, while peaks were detected on March 
29 and April 1 1 . Near-IR data in the JHK bands were taken at 
Campo Imperatore and Roque de los Muchachos (Liverpool). 
Millimetric flux densities at 345 and 230 GHz came from the 
Submillimeter Array (SMA) on Mauna Kea. Centimetric radio 
data were acquired at Medicina (22 and 8 GHz), Noto (43 GHz), 
and UMRAO (14.5, 8.0, and 4.8 GHz). In Fig. 1 the source ra- 
dio behaviour in different bands is compared to the optical one 
(top panel). The light curve at high radio frequencies (230-345 
GHz) suggests that the mechanism producing the flaring events 
observed in the optical band in the second half of February and 
in late March-April 2008 also interested the millimetric emit- 
ting zone, with some delay. An estimate of this delay is ham- 
pered by the limited data sampling. At lower radio frequencies 
(22-43 GHz) a hint of flux increase is visible in the second part 
of the light curve, while the radio flux at 5-15 GHz shows no 
trend. This suggests that the jet regions that are responsible for 
the emission at the longest radio wavelengths are not affected by 
the flaring mechanism. 



2 Calar Alto data were acquired as part of the MAPCAT (Monitoring 
AGN with Polarimetry at the Calar Alto Telescopes) project. 
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Fig. 2. Optical light curve of PKS 1510-089 (top 
panel) obtained by the GASP-WEBT during the 
period January-April 2008 compared to its radio 
flux densities at different frequencies. The vertical 
bar indicates the time of the Swift observations. 
The yellow shaded region marks the period also 
covered by the AGILE observation. 



5. Discussion and conclusions 

In the last two years, PKS 1510-089 showed a very high ac- 
tivity in y-ray band with several flaring episodes (see Pucella et 
al. 2008, Tramacere 2008, Ciprini et al. 2009, D'Ammando et 
al. 2009, Pucella et al. 2009, Vercellone et al. 2009 and Cutini 
et al. 2009). During the period January-April 2008, the source 
showed a high variability over all the electromagnetic spectrum 
from radio to y-rays, with several flaring episodes in the optical 
band and a rapid and intense flare detected in the y-ray band in 
mid-March. 

5. 1 . Modelling the Spectral Energy Distribution 

Figure [3] shows the spectral energy distribution for the AGILE- 
GRID observation period 2008 March 17-21, including quasi- 
simultaneous optical and radio data by GASP, and UVOT and 
X-ray data by S wift. Since the source brightness over the three 
days of UVOT observations remained stable, we built a unique 
SED for the whole period, including contemporaneous data at 
other frequencies, in particular radio-to-optical data from the 
GASP. The optical and near-IR data were acquired exactly in 
the period of the UVOT observations: one 7?-band datum from 



Roque (KVA) and J, H, and K data from Campo Imperatore. 
The UVOT and GASP magnitudes were corrected for Galactic 
extinction by adopting Ag = 0.416 mag, and deriving the values 
in the other bands according to Cardelli et al. (1989). To con- 
vert magnitudes into fluxes, we assumed the zero-mag fluxes by 
Poole et al. (2008) and Bessel et al. (1998). 

Radio light curves are less well-sampled than the /?-band one 
(see Fig. 2); the high-frequency radio data (43-345 GHz) shown 
in the SED were taken within a week from the UVOT obser- 
vations, while the low-frequency data points were obtained by 
interpolating between the closest data preceeding and follow- 
ing the UVOT observations. This is justified by the smooth be- 
haviour of the low-frequency radio light curves. 

The dip in the SED corresponding to the UVW1 frequency 
must be regarded with cautious, since it is also found for 
other blazars with different redshift and could be systematic. 
Observations performed by S wift and the GASP in March and 
June 2007, when the source was nearly at the same bright- 
ness level, showed the same shape in the near-IR-to-UV part of 
the SED. A similar trend characterized the optical-UV SED of 
August 2006 shown by Kataoka et al. (2008). We notice that 
the shape of the SED in the optical band may be affected by the 
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Fig. 3. Spectral Energy Distribution of PKS 1510-089 for the AGILE-GRID observation of 17-21 March 2008, including quasi- 
simultaneous GASP radio-to-optical data, the S vw'/f/UVOT data of 20-22 March and the S wiftfiW data of 20 March. The dotted, 
dashed, dot-dashed and double-dot-dashed lines represent the accretion disk black body, the SSC, the ECD and the ECC radiation, 
respectively. 



flux contribution of broad emission lines, including the little blue 
bump (Neugebauer et al. 1979, Smith et al. 1988). 

For the SED we used the S wift/XRT data collected on 20 
March, the observation closest to the y-ray flare and during 
which the higher X-ray flux was observed. 

In order to model the spectral energy distribution we 
used a homogeneous one-zone synchrotron self Compton 
(SSC; Marscher & Gear 1985, Maraschi et al. 1992, Bloom 
& Marscher 1996) model, plus the contribution of external 
Compton scattering of both direct disk radiation (ECD; Dermer 
et al. 1992) and photons from the broad line region (BLR) clouds 
(ECC; Sikora et al. 1994). The strong thermal features usu- 
ally observed in FSRQs (and in this blazars in particular, see 
Neugebauer et al. 1978, Smith et al. 1988) at optical/UV fre- 
quencies suggest that the environment is rich of soft photons 
produced by the accretion disk and/or reprocessed by the BLR. 
This implies that the energy density of the external soft radiation 
is much higher than that of the synchrotron radiation, therefore 
during the y-ray flares in FSRQs the most important processes 
are the ECC and ECD and the y-ray photon index could be de- 
termined by the dominant contribution of the two. 

We consider a moving spherical blob of radius R, filled by 
relativistic electrons and embedded in a random magnetic field. 
We assume that the electron energy density distribution is de- 
scribed by a broken power law: 



(y/7b) Pl + (yhb) Pl 

where y is the electron Lorentz factor assumed to vary between 
j m i n and y max , p\ and pi are pre- and post-break electron distri- 
bution spectral indexes, respectively, and jb is the break energy 
Lorentz factor. We assume that the blob contains a comoving 
random average magnetic field B with a bulk Lorentz factor F 



Parameter 


Value 


Units 


Pi 


2.2 




Pi 


4.6 




Tmin 


30 




7b 


290 






5200 




K 


75 


cirr 3 


R 


10 


10 15 cm 


B 


3.5 


G 


6 


20.26 




u 


5 


10 45 ergs-' 





2.86 


degrees 


r 


18 





Table 2. Parameters for the model used to explain the SED of 
PKS 1510-089 during the y-ray flare of March 18-19, 2008. 



at an angle with respect to the line of sight. The relativistic 
Doppler beaming factor is then 6 = \T{\-f3cosff)Y l , and K is the 
normalization density parameter into the blob. 

We have chosen an angle of view of 0.05 rad in agreement 
with both the apparent jet velocities derived from multiepoch 
Very Long Baseline Array (VLBA) observations of the source 
(Homan et al. 2001 ; Wardle et al. 2005; Jorstad et al. 2005; Lister 
et al. 2009) and the value used by Kataoka et al. (2008). 

The short time variability observed in y-ray band constrains 
the size of the emitting region to R < cAt var d /(1+z) = 3.86 x 
10 16 cm, where Af rar is the observed variation time. An accretion 
disk characterized by a black body spectrum with a luminosity 
of 5 x 10 45 erg s _1 , as estimated with UV observations by Pian 
& Treves (1993), at 0.05 pc from the blob is assumed as one of 
the sources of external target photons. We also assumed a BLR 
at 0.2 pc, reprocessing 10% of the irradiating continuum. The IC 
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spectra derived from the approximation of the BLR radiation as a 
black body reproduces quite well more refined spectra calculated 
by Tavecchio et al. (2008) taking into account a more accurate 
shape of the BLR. 

Assuming a model with synchrotron, SSC and EC compo- 
nents plus the contribution from the accretion disk radiation, the 
spectral energy distribution of mid March 2008 can be well rep- 
resented with input parameters summarized in Table 2. In the 
choice of the parameters we were guided by the knowledge of 
the angle of view, the disc luminosity and the simultaneous ob- 
servations of the synchrotron and IC peak regions (therefore of 
the synchrotron and IC peak frequencies and luminosities). In 
addition, the minimum variability timescale gives an indication 
of the size of the emitting region. However, even if these quan- 
tities are quite well tightly contrained by fitting the whole SED, 
the choice of some parameters is not unique because the con- 
temporaneous presence of the synchrotron, SSC, ECC and ECD 
components leads to a possible partial degeneration of the pa- 
rameters. 

5.2. X-ray spectral evolution 

The spectral evolution detected in X-rays by Swift in just two 
days, soon after the y-ray flaring episode is another hint of the 
rapid change of activity of this source. Usually in FSRQs as PKS 
1510-089 only little variability is observed in X-rays on short 
timescales from hours to days, and also on longer timescales 
the X-ray spectral shape is almost constant with only little vari- 
ations. The photon indices measured by Sw;/?/XRT in March 
2008, in particular during the first observation, tend to be lower 
than usually observed in FSRQs (<F[o.i-2.4]> = 1 .76 + 0.06 and 
<r[2~io]> = 1-65 + 0.04, Donato et al. 2001) and are more simi- 
lar to those observed in some high redshift quasars (such as RBS 
315, Tavecchio et al. 2007, and Swift J0746+2548, Sambruna 
et al. 2007). The hard photon indices of high redshift blazars 
could be interepreted in terms of absorption by warm plasma in 
the region surrounding the source, in agreement with a scenario 
where in the early evolution phases the quasars are substantially 
obscured by gas subsequently expelled from the host galaxy by 
powerful winds (see e.g. Fabian 1999). However, considering the 
low redshift, this interpretation is unlikely for this source. 

Instead the X-ray spectral evolution observed by S wift/XRT 
could be due to the contamination of an additional component 
below ~2 keV: the soft X-ray excess. In fact previous observa- 
tions with Chandra (Gambill et al. 2004) and Suzaku (Kataoka 
et al. 2008) seem to indicate the presence of the soft X-ray ex- 
cess in the spectrum of PKS 1510-089. The soft X-ray excess 
is an emission in excess of the extrapolation of the power law 
component dominating at higher energies, but the origin of this 
excess in AGNs is still an open issue (see e.g. D'Ammando et 
al. 2008b for a detailed discussion). In the past, it was often 
associated with the thermal emission of the accretion disc and 
then related to the big blue bump. However, it has recently been 
shown that modelling the soft X-ray excess in non-blazar AGNs 
with a thermal component yields a disc temperature remarkably 
constant, around 0.1-0.2 keV, regardless the central black hole 
mass and luminosity (Gierlinsky & Done 2004, Crummy et al. 
2006). Also in Kataoka et al. (2008) the soft X-ray excess is ten- 
tatively described by a black body with temperature kT ^ 0.2 
keV. This result is difficult to explain in any model for the soft 
excess related to disc continuum emission, as in any disc model 
the temperature is expected to vary with both the black hole mass 
and the accretion rates. 



1.2 



1.0 - 




17.5 17.0 16.5 16.0 15.5 15.0 

R 

Fig. 4. B—R color index versus 7?-band magnitude for PKS 1510- 
089 obtained with archive data of the Torino Observatory. 



For FSRQs one possible theoretical explanation is that the 
soft X-ray excess is a bulk Comptonization feature produced by 
cold plasma accelerated in a jet (Celotti et al. 2007), even if until 
now this feature has never been positively observed; in BL Lac 
objects instead the radiative environment is too weak to produce 
the soft X-ray excess via bulk Compton and the soft X-ray ex- 
cess is likely related to the high energy tail of the synchrotron 
emission. 

The change of photon index observed during the S wift/XRT 
observations could be due to the fact that the spectral shape of the 
inverse Compton component in X-ray remains roughly constant, 
but the amount of contamination from the soft excess emission 
varies. The contribution of the soft X-ray excess would be more 
important when the source gets fainter, affecting the spectrum at 
higher energies. Unfortunately the brief exposure of the Swift 
observation does not allow a detailed spectral modelling of this 
feature. 

A possibility for the origin of this hard power law in PKS 
1510-089 is that the photon index observed in X-rays is due 
to the combination of the synchrotron self Compton and exter- 
nal Compton emission and therefore due to the mismatch of the 
spectral slopes of these two components, not to the presence of a 
real soft X-ray excess. This is the solution that the data presented 
in this paper would favour. In this context, the spectral evolution 
during the three S wift/XRT observations could be due to the 
change of contribution of one of the two components and there- 
fore to a different variability of the SSC and EC components. 

5.3. Thermal emission components 

Even if the SED of the blazars are usually dominated by the 
beamed non-thermal jet radiation, some of them show the signa- 
ture of features Seyfert-like as the little blue bump and the big 
blue bump. The little blue bump is usually observed in quasars 
between ~ 2000 and ~ 4000 A in the rest frame and it is likely 
due to the contribution of Fell and Mgll emission lines and the 
Balmer continuum produced in the Broad Line Region (Wills et 
al. 1985). The big blue bump instead is associated with a rise 
in the UV band commonly interpreted as thermal emission from 
the accretion disc (see e.g. Laor 1990). Evidences of these ther- 
mal components have been found in other quasar-like blazar as 
3C 273 (Grandi and Palumbo 2004; Tiirler et al. 2006), 3C 279 
(Pian et al. 1999) and 3C 454.3 (Raiteri et al. 2007). 
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The presence of the emission by the accretion disk is con- 
sistent with the scenario in which the seed photons for the IC 
producing the y-ray s are external to the relativistic jet but usu- 
ally it is not observed because hidden by the beamed variable 
synchrotron emission. The fact that the synchrotron component 
of PKS 1510-089 peaks around 10 13 Hz (see Bach et al. 2007, 
Nieppola et al. 2008) allows us to observe these thermal features 
in this source. In fact around 10 15 Hz a rising emission is vis- 
ible in the spectrum and it is likely a manifestation of the big 
blue bump produced by the accretion disk, as already discussed 
for this source by Malkan & Moore (1986) and Pian & Treves 
(1993); moreover a hint of the presence of the little blue bump 
seems to appear in the SED of the source at 10 14 5 Hz. 

Given the redshift of PKS 1510-089, the Ha, H/3, Fell and 
Mgll lines mostly contribute to the observed spectrum between 
10 14 - 2 and 10 14 8 Hz, and together with the disc emission, could 
explain the excess of emission observed around 10 14 5 Hz and 
not modelled from the other components represented in the SED. 
Moreover the presence of these non-jet components in the blue 
part of the spectrum of this blazar has already been observed 
by Neugebauer et al. (1979) and Smith et al. (1988) and it is in 
agreement with the redder-when-brighter behaviour shown by 
the B - R index versus /?-band plot of Fig. 4. The plot has been 
obtained with archive data stored at the Torino Observatory. 

5.4. Energetics and alternative model 

Finally to estimate the energetics of PKS 1510-089 we compute 
the isotropic luminosity in the y-ray band, comparing with the 
Eddington and bolometric luminosity and the total power trans- 
ported by the jet. For a given source with redshift z, the isotropic 
emitted luminosity in the energy band e is defined as: 



4nFd 2 (z) 
(1 +z) (1 - a ? ) ' 



(2) 



where, in our case, e is the y-ray energy band with L m ; n = 
lOOMeV and E m . dx = lOGeV, a y = T - 1, F is the y-ray en- 
ergy flux between L m ; n and L max calculated from the photon flux 
F y (E > 100 MeV) as suggested by Ghisellini et al (2009): 



F = 1.6 x 10- 12 -^^[100 1- ^ 
1 — a y 



1] 



The luminosity distance is given by 
diizuzi) = (1 +Z2? x r[E(z)]- l dz, 

1 + Z2 J Zl 

where Zi = 0, Z2 = z src and 

E(z) = VOm(l + z) 3 + (1 - fi m - Q A )(1 + z) 2 + Q A , 



(3) 



(4) 



(5) 



where Ho is the Hubble constant, O m and Qa are the contribution 
of the matter and of the cosmological constant, respectively, to 
the density parameter. Using a luminosity distance d L = 1915 
Mpc and the average y-ray flux observed by the AGILE-GRID 
during 17-21 March 2008, we obtain for PKS 1510-089 (z = 
0.361) an isotropic luminosity Ly° = 5.3 x 10 47 erg s _1 . 

The power carried by the jet in the form of magnetic field 
(L B ), cold protons (L p ), relativistic electrons (L e ) and produced 
radiation (L ra ^), are: 



L p = 7i R 2 F 2 c j [N(y) m p c 2 dy] = 3.6 x 10 45 erg s" 1 



nR 2 Y 2 c j [N(y)ym e c 2 dy] = 1.5 x 10 44 ergs-> (7) 



L B =nR 2 T 2 cU B = 1.5 x 10 45 erg s" 1 
Lr ad - L iso F 2 /6 4 = 5.3 x 10 45 erg s" 1 



(8) 
(9) 



where Ub is the magnetic energy density. Therefore the total 
power transported by the jet is P = Lb + L p + L e + L ra d = 1.1 x 
10 46 ergs _1 . 

Assuming for the source a black hole mass Mbh = 4.5 x 10 8 
M Q (Woo & Urry 2002), we obtain an Eddington luminosity 



Lva& = *JI^Lm bh = 5.7 xlO 46 ergs- 



(10) 



(6) 



to be compared with the bolometric luminosity L/, i = 2.4 x 
10 46 erg s" 1 reported in Woo & Urry (2002). 

An alternative theoretical model was recently proposed by 
Kataoka et al. (2008) to interprete the data of PKS 1510-089 
collected during August 2006, with the high energy emission 
originated by the Comptonization of infrared radiation produced 
by the molecular torus surrounding the central engine and sug- 
gesting that the soft X-ray excess could be produced by the 
IC scattering of external photons by a population of cold elec- 
trons, as discussed by Begelman et al. (1997) and Celotti et al. 
(2007). An accurate theoretical interpretation of the SED is be- 
yond the scope of this paper; however, our data seem not to 
rule out this alternative model. Infact, we have not simultane- 
ous observations in FIR band that can confirm the excess de- 
tected by IRAS (Tanner et al. 1996), interpreted by Kataoka et 
al. as due to dust radiation from the nuclear torus and assumed as 
main source of seed photons for the IC mechanism. Moreover, 
the bulk Comptonization feature should not be easily observable 
during a high activity state of the source, such as that observed 
in mid-March 2008, because overwhelmed by the SSC and ECD 
emission and therefore this is not the best situation to test this 
hypothesis. 

Further X-ray observations with XMM-Newton and S uzaku, 
simultaneously with the optical monitoring by means of REM 
Telescope and WEBT Consortium, could give important indi- 
cation of the emission mechanisms involved in this source, in 
particular of the real nature of the soft X-ray excess, the pres- 
ence in the spectrum of Seyfert-like features and the existence 
of the bulk Comptonization feature or not. 

Finally, with two y-ray satellites, AGILE and Fermi, in orbit 
at the same time we will be able to study in detail the source 
behaviour at high energies on long time scale, even if a wide 
multiwavelength coverage is necessary to achieve a complete 
understanding of the structure of the jet, the origin of the seed 
photons for the inverse Compton process and all the emission 
mechanisms working in this blazar. 
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